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ABSTRACT 



yA/ms. The aim is to explore the interstellar medium around the dust bubble N131 and search for signatures of star formation. 
Methods. We have performed a multi wavelength study around the N131 with data taken from large-scale surveys of infrared obser- 
vation with online archive. We present new observations of three CO 7=1-0 isotope variants from Purple Mountain Observatory 
13.7 m telescope. We analyzed the distribution of the molecular gas and dust in the environment of the N131. We used color-color 
diagrams to search for young stellar objects and to identify ionizing star candidates. 

Results. The kinematic distance of ~8.6 kpc has been adopted as the distance of the bubble N131 from the Sun in this work. We 
find a ring of clouds in CO emission coincident with the shell of N131 seen in the Spitzer telescope images, and two giant elongated 
molecular clouds of CO emission appearing on opposite sides of the ringlike shell of N131. There is a cavity within the bubble at 1.4 
GHz and 24 fim. Seven IRAS point sources are distributed along the ringlike shell of the bubble N131. Fifteen ionizing stars and 63 
YSO candidates have been found. The clustered class I and II YSOs are distributed along the elongated clouds in the line of sight. 

Key words, infrared: stars — stars: formation — ISM: bubbles — HII regions — radio lines: ISM 



1. Introduction 



IChurchwell et all (l2006i l2007l) have detected and cataloged 
about 600 midinfrared dust (MIR) bubbles between longitudes - 
60° and +60°. The IR dust bubbles may be produced by ionizing 
O- and/or B-type stars, which are located inside the bubble. The 
ultraviolet (UV) radiation from ionizing stars may heat dust and 
ionize the gas to form an ex panding bubble shell (Wa tson et al.l 
2008I). [Simpson et all (120121) present a new catalog of 5106 in- 
frared bubbles created through visual classification via the on- 
line citizen science website "The Milky Way Project' Q, which 
provides a crowd-sourced map of bubbles and arcs in the Milky 
Way, and will e nable better statistical analy sis of Galactic star- 
forming sites. iBeaumont & Williamsl (1201 Oh report CO J = 3 - 
2 maps of 43 Spitzer identified bubbles, which suggests that ex- 
panding shock fronts are poorly bound by molecul ar gas and that 
cloud compression by these shocks may be limited. IWatson et a"D 
(2008) present an analysis of wind-blown, parsec-sized, midin- 
frared bubbles, and associated star formation, and suggest that 
more than a quarter of the bubbles may have triggered the for- 
mation of massive objects. 

A few individ ual bubb les have been stud i ed well, such as 
N22 dJi et al.ll2012l). N49 (IWatson et al.ll2008l: iDeharveng et al.l 
120101; IZavagno et al.1 l2010h . N6 5 dPetriella et al l l2010h . N68 
(IZhang & Wangl2012bl) . and S51 (IZhang & Wangl2012al) . There 
are many models and observations to e xplain the dusty wind- 
blown bubbles, such as bubble N49 of lEverett & Churchwelll 
(120101) . Recently, we have reported an expanding ringlike shell 
of the bubble S51, which shows a rare front side located within 
the shell in the line of sight, by employing 13 CO and C 18 J 



= 1-0 emission lines of the Mopra Telescope (Zhang & Wang 
I20l2^ . We also investigated the star formation around the bub- 
ble N68, which suggests th at the massive star form ation at the 
ringlike shell is very active (Zhang &Wang 2012b). 

Toward the bubble N131, in this work, we carried out new 
observations of the / = 1-0 transitions of 12 CO, 13 CO, and 
C 18 using the telescope of the Purple Mountain Observatory 
(PMO) at Qinghai province in China, we mainly report several 
molecular clumps associated well with the IR dust bubble N13 1 . 
We also aim to explore its surrounding ISM and search for star 
formation spots. We describe the data in Sect. |2 the results and 
discussion are presented in Sect. [3J Sect. [4] summarizes the re- 
sults. 



2. Observation and data processing 

2.1. The online archive 

The data of t he online archive used in this work in - 
clude GLI MPSE ([Be njamin e t al.l2003l : lchurchwell et al.l2009h . 
MIPSGAL dCarev et al.ll2009h. the Two Micr on All Sky Survey 
(2MA SSfl (SkrutskieeLaLl l20Q£|), IR AS dNeugebauer et al.l 
I1984 . and NVSS dCondon et al.lll998l) . GLIMPSE is an MIR 
survey of the inner Galaxy performed with the Spitzer Space 
Telescope. We used the mosaicked images from GLIMPSE and 
the GLIMPSE Point-Source Catalog (GPSC) in the Spitzer- 
IRAC (3.6, 4.5, 5.8, and 8.0 u m). IRAC has an angular reso- 
lution between 1.5" and 1.9" (Faz io et al.l 12004: We rner et all 



http ://w w w. milky wayproj ect . org 



2 2MASS is a joint project of the University of Massachusetts 
and the Infrared Processing and Analysis Center/California Institute 
of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 
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Table 1. Parameters of the nine spectra 



Clump 


R.A.(J2000) 


DEC.(J2000) 


V 12 C 


V 13 CO 


AVi2 co 


AVi3 co 


T i2 C o 


T 13 C 




h m s 




(kms- 1 ) 


(km s~ 2 ) 


(kms -1 ) 


(km s" 2 ) 


(K) 


(K) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


A 


19 51 55.681 


26 21 29.87 


-8.91±0.01 


-9.25+0.03 


2.63+0.02 


2.14+0.06 


11.40+0.25 


2.04+0.13 


B 


19 52 45.568 


26 25 54.58 


-11.89+0.02 


-12.15+0.03 


2.65+0.05 


2.05+0.07 


7.19+0.31 


2.38+0.16 


C 


19 53 01.891 


26 28 06.93 


-12.24+0.04 


-12.53+0.08 


2.62+0.09 


2.53+0.17 


4.13+0.26 


0.95+0.16 


D 


19 51 32.383 


26 20 22.06 


-9.93+0.02 


-10.15+0.10 


2.30+0.05 


2.23+0.29 


6.07+0.34 


0.95+0.19 


E 


19 52 18.309 


26 27 02.22 


-8.09+0.03 


-8.29+0.05 


1.55+0.09 


1.63+0.17 


5.22+0.23 


1.02+0.16 


F 


19 52 51.126 


26 19 23.40 


-9.80+0.03 


-9.84+0.12 


1.77+0.06 


2.05+0.37 


5.25+0.36 


0.62+0.16 


G 


19 52 35.316 


26 18 01.05 


-10.64+0.01 


-10.68+0.02 


1.97+0.03 


1.33+0.05 


7.35+0.31 


2.24+0.12 


H 


19 52 16.594 


26 15 45.75 


-12.15+0.03 


-12.36+0.10 


1.82+0.08 


1.21+0.22 


2.96+0.27 


0.50+0.11 


I 


19 52 02.983 


26 16 21.05 


-8.31+0.04 


-8.33+0.09 


3.85+0.10 


2.72+0.28 


3.15+0.29 


0.73+0.15 



Table 2. Derived parameters of the CO molecular clumps 



Clump 


Area 


Int.Ti2 C0 


Int.Ti3 co 


T ex (max) 


T ex (mean) 


N H ,(max) 
(10 21 cm" 2 ) 


Nh 2 (mean) 
(10 21 cm" 2 ) 


Mh 2 


n H 2 




(arcmin 2 ) 


(arcmin 2 K) 


(arcmin 2 K) 


(K) 


(K) 


(Mo) 


(cm" 3 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


AD 


36.59 


554.34 


102.31 


29.15 


11.54 


15.39 


3.36 


1438 


72 


BC 


32.23 


491.05 


122.32 


23.17 


12.42 


21.10 


4.71 


1775 


107 


A 


13.72 


305.22 


50.05 


29.15 


14.56 


15.39 


5.07 


814 


177 


B 


12.63 


233.59 


58.97 


23.17 


14.73 


21.10 


5.69 


841 


207 


C 


6.75 


118.92 


30.93 


21.94 


17.12 


17.12 


8.48 


669 


423 


D 


2.18 


37.48 


6.66 


20.90 


15.12 


8.90 


5.00 


127 


439 


E 


3.70 


31.23 


7.71 


17.10 


9.28 


3.71 


1.25 


54 


84 


F 


3.05 


40.64 


7.53 


17.79 


12.90 


4.88 


2.12 


76 


157 


G 


3.70 


50.50 


13.38 


18.95 


11.19 


8.42 


2.53 


109 


170 


H 


1.31 


10.01 


2.41 


11.89 


8.28 


1.75 


0.58 


9 


66 


I 


2.18 


27.49 


5.15 


13.08 


8.03 


3.35 


1.06 


27 


93 



Table 3. IRAS point sources around the bubble N131 



Name 


IRAS Source 


R.A.(J2000) 


DEC.(J2000) 


Fu 


^25 


^60 


^100 


Lir 


T d 






h m s 




Jy 


Jy 


Jy 


Jy 


Lo 


K 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


IR1 


19498+2621 


19 51 58.18 


26 29 46.36 


0.31 


0.43 


3.14 


33.82 


1855 


19.47 


IR2 


19499+2613 


19 52 01.47 


26 21 09.56 


5.54 


22.02 


216.90 


322.80 


40873 


32.52 


IR3 


19500+2617 


19 52 09.90 


26 24 56.10 


0.45 


0.27 


3.23 


322.80 


13199 


13.41 


IR4 


19501+2607 


19 52 12.62 


26 15 03.28 


0.44 


1.12 


9.84 


31.69 


2627 


25.78 


IR5 


19502+2606 


19 52 19.94 


26 14 19.75 


4.14 


1.70 


2.99 


31.51 


3717 


19.56 


IR6 


19502+2618 


19 52 21.28 


26 26 41.84 


0.28 


0.27 


2.55 


53.40 


2520 


17.16 


IR7 


19503+2617 


19 52 26.01 


26 25 28.14 


5.96 


2.60 


2.72 


43.68 


5149 


18.01 


IR8 


19505+2610 


19 52 40.28 


26 18 04.06 


0.30 


0.22 


3.38 


23.33 


1429 


21.40 



12004 . MIPSGAL is a survey of the same region as GLIMPSE, 
using the MIPS instrument (24 and 70 yum) on Spitzer. The 
MIPSGAL resolution is 6" at 24 jim. The IRAS Point Source 
Catalog consists of 245889 sources found and verified by the 
IRAS (InfraRed Astronomy Satellite) at 12, 25, 60, and 100 
yum. The NRAO VLA Sky Survey (NVSS) is a 1.4 GHz con- 
tinuum survey covering the entire sky north of -40° declination 
(ICondon et alll998h : the NVSS survey has a noise of about 0.45 
mJy beam -1 . 

2.2. The CO data of Purple Mountain Observatory 

Our CO observations were made during May 2012 using the 
13.7-m millimeter telescope of Qinghai Station at the Purple 
Mountain Observatory at Delinghaj. We used the nine-pixel ar- 
ray receiver separated by -180". The receiver was operated in 
the sideband separation of single sideband mode, which allows 
for simultaneous observations of three CO J = 1 - isotope 
variants, with 12 CO in the upper sideband (USB) and 13 CO and 
C 18 O in the lower sideband (LSB). The half-power beam width 
(HPBW) is 52"±3", and the main beam efficiency is -50% at 
-110 GHz. The pointing and tracking accuracies are better than 
5". The typical system temperature during our runs was around 
110 K and varies by about 10% for each beam. A fast Fourier 



3 http://www.dlh.pmo.cas.cn/ 



transform (FFT) spectrometer was used as the back end with a 
total bandwidth of 1 GHz and 16384 channels. The velocity res- 
olution is about 0.16 km s" 1 at ~1 10 GHz. 

On-the-fly (OTF) observing mode was applied for mapping 
observations. The antenna continuously scanned a region of 
20'x20' centered on R.A.(J2000) = 19*52 m 09 5 .61, DEC.(J2000) 
= 26°22 , 13 ,, .4 with a scan speed of 20" s _1 . The OFF posi- 
tion was chosen at R.A.(J2000) = 19M7 m 24*.00, DEC.(J2000) = 
28°29 / 24".0, where there is extremely weak CO emission base d 
on the CO survey of the Milky Way ( Dame et al . 1987l l200Tb . 
The rms noise level was 0.2 K in main beam antenna tempera- 
ture T^ for 12 CO (1-0), and 0.1 K for 13 CO (1-0) and C 18 (1-0). 
The OTF data were then converted to 3-D cube data with a grid 
spacing of 30". The IRAM software package GILDAS0 and the 
software package MIRIAE0 were used for the data reduction. 



3. Analysis and the results 

3.1. The dimension and distance of the bubble N131 

We selected the IR d ust bubble N131 from the catalog of 
IChurchwell et al.l (l2006h . They suggest that the Nl 31 is a com- 
plete (closed ring) IR dust bubble centered on /* = 63.084, b* 



4 http://iram.fr/IRAMFR/GILDAS/ 

5 http ://w ww. cf a. harvard, edu/sma/miriad/ 
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Table 4. Ionizing star candidates within bubble 



GLIMPSE Desig. 


m.e 


m A .5 


™5.8 


^8.0 


mj 


m H 


m Ks 


Ay 


Mj 


M H 


Mk s 


O-type star 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


G063. 0696-00. 3941 


11.34 


11.39 


11.29 


11.23 


12.63 


11.89 


11.60 


7.06 


-4.04 


-4.02 


-3.86 


07V-06.5V 


G063.08 19-00.3947 


11.01 


10.94 


10.88 


10.89 


12.65 


11.77 


11.29 


9.26 


-4.64 


-4.52 


-4.42 


05.5V-05V 


G063. 09 18-00. 3962 


11.75 


11.62 


11.43 


11.58 


13.60 


12.59 


12.08 


10.09 


-3.91 


-3.85 


-3.72 


07.5V 


G063.0958-00.4059 


11.75 


11.82 


11.65 


11.67 


13.70 


12.51 


12.08 


10.29 


-3.88 


-3.97 


-3.75 


07.5V-07V 


G063.0962-00.3644 


11.59 


11.54 


11.51 


11.42 


13.04 


12.07 


11.70 


8.76 


-4.11 


-4.14 


-3.95 


07V-06.5V 


G063. 1020-00.4375 


11.55 


11.60 


11.48 


11.43 


13.51 


12.32 


11.87 


10.45 


-4.12 


-4.18 


-3.98 


07V-06V 


G063. 1050-00.4143 


12.16 


12.16 


11.96 


12.26 


14.20 


12.95 


12.53 


10.51 


-3.43 


-3.57 


-3.32 


09V-08.5V 


G063. 1080-00.4197 


10.82 


10.90 


10.71 


10.71 


12.61 


11.49 


11.05 


10.01 


-4.89 


-4.93 


-4.74 


04V-04V 


G063.ll 18-00.4060 


10.73 


10.65 


10.60 


10.46 


12.31 


11.37 


10.94 


9.05 


-4.92 


-4.89 


-4.74 


04V-04V 


G063. 1127-00.4278 


12.18 


12.16 


11.86 


12.00 


14.11 


12.89 


12.45 


10.51 


-3.53 


-3.62 


-3.40 


09V-08.5V 


G063. 1138-00.4385 


12.23 


12.13 


12.08 


11.99 


14.02 


12.95 


12.54 


9.58 


-3.35 


-3.40 


-3.21 


09.5V-09V 


G063. 1143-00.3893 


11.09 


11.13 


11.02 


11.04 


12.43 


11.57 


11.26 


7.81 


-4.44 


-4.47 


-4.29 


05.5V-05V 


G063. 1168-00.4392 


10.70 


10.79 


10.61 


10.62 


12.57 


11.44 


11.02 


9.94 


-4.90 


-4.97 


-4.77 


O4V-03V 


G063. 1187-00.4133 


12.00 


11.97 


11.73 


11.93 


13.91 


12.71 


12.25 


10.55 


-3.74 


-3.81 


-3.60 


O8V-07.5V 


G063. 125 1-00.4258 


10.71 


10.73 


10.54 


10.57 


12.56 


11.43 


11.00 


10.03 


-4.94 


-5.00 


-4.80 


O4V-03V 




-20 -15 -10 -5 -20 -15 -10 -5 -20 -15 -10 -5 
Velocity (km/s) Velocity (km/s) Velocity (km/s) 



Fig. 1. 12 CO (red line), 13 CO (black line), and C 18 (green line) spectra at the peaks of the molecular clumps from A to I. The 
brightness temperature of each 13 CO spectrum is multiplied by 2. 



= -0.395, with an inner short radius r* n = 5.46', an inner long 
radius R* n = 6.18', and an eccentricity of the ellipse e^ l3l = 
0.47. By comparing the IR emission with the integrated inten- 
sity of CO emission, the dimensions about the ring of cloud 
are respectively r in = 5.20', R iW = 6.00', and ^Ni3i = 0-50 cen- 
tered on / = 63.095, b = -0.404, or R.A.(J2000) = 19 /z 52 m 2P.5, 
DEC.(J2000) = +26°21 , 24 ,, .0. Here, r in and R in are the semimi- 
nor and semimajor axes of the inner ellipse, respectively. By an- 
alyzing H2CO absorption (Vh 2 co = 22.6+0.1 km s" 1 ) against 
the UC HIT region contin uum emission (Vmioa = -9.3+2.3 km 
s" 1 ), IWatson et al.l (120031) resolved the distance ambiguity to- 
ward G63.05-0.34, which lies on the "far" kinematic distance 
8.6+1 kp c - In fact, there is no distance ambiguity toward N131, 
because the derived "near" kinematic distance is negative. The 
G63.05-0.34 is located at the position of the IRAS 19499+2613, 
which is correlated with CO molecular cl ump A (Vco = — 10-5 
km s" 1 ) in this work. lWatson et al. (2010 ), however, ado pted the 
velocity of H 2 CO (-22.6 km s _1 ) in lWatson et"afl (l2003h to ob- 
tain a kinematic distance of ~2.4 kpc toward the bubble N131. 
Actually, there is another CO velocity component at -25.0 km 
s" 1 , which is consistent with H2CO velocity at -22.6 km s" 1 , but 



not correlated with the ringlike shell of this bubble. This velocity 
component of ~ 22. 6 km s" 1 may belong to the foreground of the 
bubble N131. Therefore, we adopt the kinematic distance Dni3i 
= 8.6 kpc as the distance of the bubble N131. The dimensions of 
the inner short radius and the inner long radius are D r . n = 13.0 pc 
and Dr. w = 15.0 pc, respectively. 



3.2. Parameters of the spectra and molecular clumps 

Figure [T] shows several spectra for 12 CO, 13 CO, and C 18 at the 
peak locations of nine molecular clumps from A to I (indicated 
in Fig. [2]) defined on the basis of 13 CO contours. For the nine 
molecular clumps, we detected strong 12 CO and 13 CO emission 
lines. At any position toward the N131, however, we did not de- 
tect any C 18 emission signal of more than 3<x. We reported 
spectral information of the nine molecular clumps in Table [T] 
and some derived parameters of the molecular clumps in Table O 
In Table[TJ column (1) lists the molecular clump name; columns 
(2) - (3) list the equatorial coordinates; columns (4) - (9) list 
the LSR velocity V C o, the full width at half maximum (FWHM) 
AVco, and the peak brightness temperature Tco of each clump 
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10 15 20 25 30 35 




19 h 53 m 00 s 040 s 20 s 00 s 40 s 20 s 
R.A. (J2000) 

Fig. 2. Integrated intensity contours of the 12 CO emission clouds superimposed on the GLIMPSE 8.0 jim colorscale. The contour 
levels range from 5.35 to 32. 1 1 by 2.68 K km s" 1 . The integration range is from -14.5 to -6.5 km s" 1 . The letters from A to I indicate 
the positions of nine molecular clumps, and the area of each clump is indicated with dashed ellipse and polygon. The black symbols 
"■" indicate the positions of IRAS point sources. The white ellipse indicates the position of bubble N131, and the straight lines 
indicate the position of position-velocity diagram in Fig. [5] The beam size of 12 CO emission is given in the filled red circle. The 
unit of the color bar is in MJy sr" 1 . 



using Gaussian fitting for CO emission line. In Tabled columns 
(1) - (4) show the name, the dimension, and the integration tem- 
perature of each molecular clump, where the area of each clump 
is indicated with a dashed ellipse and polygon in Fig. [2 columns 
(5) - (8) show the maximum and mean value of the ionizing tem- 
perature and the H2 column density; columns (9) - (10) show the 
mass and density of each molecular clump. Here, we assumed 
each molecular clump is under the local thermal equilibrium 
(LTE) assumption, and u sed the theory of the radiation transfer 
and m olecular excitation (iWinnewisser et al.ll979t iGarden et al.l 
Il99ll) . The excitation temperature T ex and column density Ni3 C0 
can be calculated directly, assuming 12 CO emission to be opti- 
cally thick and the beam-filling factor to be unity. The column 
densities of H2 were obtained by adopting typical abundance ra- 
tios [H 2 ]/[ 12 CO]=10 4 and [ 12 CO]/[ 13 CO]=60 in the ISM. 

3.3. The ringlike shell 

In Fig. [3 the 8.0 jim emission colorscale outlines clearly the 
distribution of molecular filament and clumpy structure toward 
the bubble N131. The 8.0 /jm emission originates mainly in 
the poly cyclic aromatic hydrocarbons (PAHs). Inside the N131, 
the 8.0 /jm emission is much weaker than that at the ringlike 
shell. It is likely that the stellar wind from O- and/or early B- 
type stars have blown the PAHs onto the ringlike shell. In Fig. 
O the colorscale is the MIPSGAL 24 /mi emission. Generally, 
there should be strong 24 /mi emission and 1.4 GHz con- 
tinuum emission inside bu bble, such as in the dust bubble s 
S51 (IZhang&Wandl2012ah and N68 (IZhang & Wand 1201 2bl) . 



However, having nearly the same distribution as 8.0 jim emis- 
sion, the 24 /jm emission is very weak inside the bubble N131. 

Comparing 12 CO emission in Fig. [2] with 13 CO emission in 
Fig. [3] we can find that several molecular clumps were connected 
together to form the ringlike shell of the Nl 3 1 . The inner edge of 
the ringlike shell has a much sharper gradient of 12 CO emission 
than the outer edge. Based on the distribution of optically thin 
13 CO in Fig. [3] we also found seven clumps A, B, E, F, and 
I along the ringlike shell. The parameters of the clumps in Table 
[2] show that the dense cores of star formation may be forming in 
clumps. Therefore, the ringlike shell is a possible birth place of 
star formation that was triggered by the bubble N131. 

3.4. Two elongated molecular clouds 

In Fig. [21 we also found there are two giant elongated molecular 
clouds (AD and BC) of CO emission appearing on opposite sides 
of the ringlike shell of N131. Morphologically, each cloud is in 
alignment, and the central axis of each cloud nearly goes through 
the center of the N131. The densest position of each cloud is 
located at the ringlike shell of the N131, and the two clouds AD 
and BC have extensive structure outwardly. On the north and 
south of the clump A located at the ringlike shell, there also is an 
expanding tendency in the western direction. This morphology 
may be caused by the stellar wind within the bubble, but the 
possibility needs to be explored further. 

Using the channel map (Fig. |4]) of 12 CO emission to inves- 
tigate the velocity component of the bubble N131, The green 
contours in the panel of Fig. [4] show that the systematic veloc- 
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Fig. 3. Integrated intensity contours of the 13 C0 emission clouds superimposed on the MIPSGAL 24 jjm colorscale. The contour 
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filled red circle. The unit of the color bar is in MJy sr" 1 . 



ity of the bubble is about -10.5+0.5 km s" 1 , which is consis- 
tent with the velocity of the ionized gas (Vmioa = -9.3+2.3 km 
s" 1 . 1 Watson et al.l (l2003h ). The velocity range of the cloud AD is 
from -9.5 to -6.5 km s" 1 indicated with red contours, while the 
velocity of the cloud BC is from -14.5 to -1 1.5 km s" 1 indicated 
with blue contours. We argue that the clouds AD and BC are, 
respectively, redshifted and blueshifted relative to the ringlike 
shell of the bubble. 

In addition, along the central axis of the clouds AD and BC 
in Fig. [2l we made a position- velocity diagram in Fig. [5] Both 
the velocity gradient and velocity dispersion gradient are obvi- 
ous from clumps A (B) to D (C). The velocity of the clump A 
on the ringlike shell is higher than that of the clump D, which is 
consist ent with the redshifted velocity distribution of the bubble 
N68 (IZhang & Wangll2012bl) . For the cloud BC, the position- 
velocity diagram shows that clumps B and C may be two inde- 
pendent components. However, it is also possible that the molec- 
ular clumps B and C may be interacting. Figure shows inte- 
grated intensity contours of the 12 CO emission of the blueshifted 
(from -14.5 to -1 1.5 km s" 1 ) and redshifted (from -9.5 to -6.5 km 
s" 1 ) clouds. There is almost no overlap region between the blue- 
and redshifted clouds. 

3.5. Lyman continuum flux 

Figure [6] shows the 1.4 GHz NVSS continuum emission con- 
tours. The region of the continuum flux value above 2<x noise 
(or > 0.9 mJy beam -1 ) was only integrated to consider as the 
reliable HII region candidates. From Fig. it can be seen that 
there is no radio continuum emission within N131, except for a 



clump projected on the plane of the sky towards the center of 
the bubble. This fact, along with that the 24 jim emission being 
very weak inside the bubble, would indicate that hot dust and 
ioniz ed gas have been evacuated by stellar wind (Wa tson et al] 
120091) . We also found that IR2 is located a t the peak of HII re- 
gion G63.049-0.348 dGregorv et al.lll996l) . suggesting that the 
IR2 (IRAS 19499+2613) may be the ionizing star of the HII re- 
gion. 

Assuming that the only radio feature in the center of N131 
is associated with the bubble, we derive the radio flux density 
of the HII region in about 0.016 Jy. The flux was estimated by 
integrating each pixel for the signals of more than 2<x noise from 
the NVSS fit image. This value should be taken with caution be- 
cause the NVSS survey has not added the flux contribution from 
large-scale structures. The number of stellar Lyman photon, ab- 
so rbed by the gas in th e HII region, follows the relation Eq. (Q]) 
in lMezger et al.l(ll974l) 



Ni 
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4.761xl0 48 -a(v,T e r 
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a) 

where a(v, T e ) is a slow ly varying function tabulated by 
Mezger & Hende rson! (Il967l) . for effective temperature of ioniz- 
ing star T e ~ 33000 K and at radio wavelengths, a(v, T e ) ~ 1. 
Finally, we obtained Lyman continuum ionizing photons flux 
logNLyc ~ 46.73 from nebula. Assuming th e ionizing sta rs be- 
long to the 09.5 star with logN Lyc ~ 47.84 dPanagiall 1973b . and 
the estimated Nl yc is a factor of 2 lower than the expected ioniz - 
ing photon flux from these stars (Beaumont & Williams 20101) . 
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Fig. 4. Integrated intensity contours of the 12 C0 emission clouds every 1 .0 km s" 1 superimposed on the GLIMPSE 8.0 jim grayscale. 
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position of bubble Nl 3 1 . 



it is suggested that there should be about 0.16 ionizing star of 
09.5 star to ionize the ISM within the bubble N131. 

We also used the Effelsberg radio continuum data from the 
Galactic plane survey at 2695 MH z (only source co mponent) to 
derive the Lyman continuum flux dFurst et al.| [l990). The beam 
size of this survey is about 4.3', rms is 20 mK in brightness tem- 
perature, and T B /S = 2.51+0.05 [K/Jy]. Within a circle of 2.60' 
(0.5xr;„) in radius centered on Nl 3 1 , the average brightness tem- 
perature is about 0.3 K, so the intensity is about 0.120 Jy. Using 
Eq. (Q}, we obtained the Lyman continuum flux logNL yc ~ 47.63. 
After considering the underestimated factor 2, we found that 
there is only about 1.24 ionizing star of 09.5 star to ionize the 
ISM. Therefore, the Effelsberg radio continuum intensity is also 
weak like the 1 .4 GHz NVSS continuum. Maybe it is not enough 
to ionize the bubble N131 for the weak continuum intensity. 

However, we found 15 reliable ionizing star candidates 
within the bubble (see Sect. l3.7b . It is likely that hot dust and ion- 
ized gas have been evacuated by a stellar wind from the clustered 
ionizing stars, and this will lead to underestimating the Lyman 
continuum ionizing photons flux from the ionizing stars. 



3.6. IRAS point sources 

Within a circle about 10' in radius centered on N131, we found 
eight IRAS point sources, indicated in Figs. |2 [3] and [6] with the 
names IR1, IR2, and IR8. These IRAS sources, except IR1, 
are distributed around the several molecular clumps in the line 



of sight. Especially, IR4, IR6, and IR8 are correlated well with 
the molecular clumps A, H, E, and G, respectively. We obtained 
some parameters of the IRAS point sources in Table [3] Columns 
(1) - (2) list the source name; columns (3) - (4) list the equa- 
torial coordinates; columns (5) - (8) list the flux of the 12, 25, 
60, and 100 /mi, respec tively; columns (9 ) - (10) list the derived 
the infrared luminosi ty (iCasoli et al.ll 1986b and dust temperature 
( Hen ning et al.lll990l) . which are expressed as 

L IR = (20.653x J F 12 +7.538x J F25+4.578x J F 60 +1.762x J F 1 oo)xZ) 2 x0.30, 

(2) 

96 

Td " (3+^)ln(100/60)-ln(F 6 o/F 10 o)' (3) 
In the equations above, D is the distance from the Sun in kpc, 
and the emissivity index of dust particles j3 is assumed to be two. 
Based on the parameters about the infrared luminosity, dust tem- 
perature in Table [3] the eight IRAS point sources are probable 
candidates for young massive stars. 

In studying of the occurrence of maser emission from 
star-forming r e gions in the very early stages of evolution, 
Palu mbo et al.l (1 19941) found the 22 GHz H 2 Q maser emis sion 
has a flux of less than 3.1 Jy, and Ivan der Walt et aD ([1995) did 
not find any 6.7 GHz methanol maser toward the IR2 (IRAS 
19499+2613). Also, in a study surveying the occurrence of the 
22 GHz H2O m aser emission fro m bright IR sources in star- 
forming regions, iPalla et aD (Il99ll) found the peak flux of H2O 
maser is less than 2.7 Jy toward the IR4 (IRAS 19501+2607) 
near clump H. 
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for the right panel. The color bar of the grayscale for the 12 CO emission is above the diagram. The letters A, B, C, and D correspond 
to the positions of the same letters as in Fig. [2] 



3. 7. Ionizing stars and YSOs 

The GLIMPSE color-color (CC) diagram [5. 8] -[8.0] versus 
[3.6]-[4.5] in Fig.|3a) shows the distribution of class I, II, and 
III stars, which are located within a circle of about 10' in radius 
centered on N131. Here we only consi dered these sourc es with 
detection in four Spitzer-IRAC bands dHora et al.l l2008). Class 

I sources are protostars with circumstellar envelopes; class II 
sources are disk dominated objects; and class III sources refer to 
young stars above the main sequence (and contracting towards 
it), but without accretion characteristi cs such as Ha emission 
(lAUen et al. 1120041 iPetriella et alJl2010h . 

The reliable ionizing stars are mainly from the class III can- 
didates and centrally distributed in a small region within the bub- 
ble, so we only considered the sources within a circle of 2.60' 
(0.5xr in ) in radius centered on N131. Furthermore to get rid of 
the background and foreground stars, we used the 2MASS CC 
diagram (H - K s ) versus (/ - H) in Fig. |3b). Considering the 
extin ction in the Galact i c pla ne as a function of distance to the 
Sun, lAmores & Lepinel (120051) present an average Galaxy ISM 
extinction Ay =0.96 mag kpc" 1 for one model. The adopted dis- 
tance of the N131 is 8.6 kpc from the Sun, so we just consid- 
ered the sources between the extinction range 10.6 ~ 6.6 mag. 
Finally, we obtained 15 reliable ionizing star candidates within 
the bubble Nl 3 1 . In Tabled we report the 1 5 ionizing candidates 
within the bubble: column (1) specifies the GLIMPSE designa- 
tion; columns (2) - (8) list the magnitude of four Spitzer-IRAC 
bands and three 2MASS JHK S bands, respectively; column (9) 
lists the ISM extinction Ay; the absolute JHK S magnitude and 
spectral type were derived in columns (10) - (13). 

YSOs are mainly distributed around the ringlike shell of the 
bubble, so we selected YSO candidates from class I and class 

II sources in Fig. |3a). We then excluded the potential sources 
belonging to the reddening main sequence and giant stars us- 
ing the 2MASS CC diagram in Fig. Eft). These (class I and II) 
sources in Fig. |3b) only include those with simultaneous de- 
tection in the four Spitzer-IRAC bands and three 2MASS JHK S 
bands. Finally, we found 29 class I stars and 34 class II stars as 



YSO candidates, which are indicated around the bubble N131 
in Fig. [3] We can see that the clustered YSOs are distributed on 
the clouds AD and BC in the line of sight, and several YSOs are 
located near the molecular clumps G and H on the ringlike shell. 
This distribution provides some evidence for star formation trig- 
gered by the bubble. Within the bubble, there are two YSO can- 
didates, which are possibly the background or foreground stars. 
We also found that there is a good correlation between the YSOs 
and MIPSGAL 24 /mi distribution. 

4. Summary 

Based on our CO emission observations of 13.7-m PMO tele- 
scope, together with other archival data including GLIMPSE, 
MIPSGAL, 2MASS, IRAS, and NVSS, we have studied the ISM 
around the IR dust bubble N131. The main results can be sum- 
marized as follows. 

1. We found that the ringlike shell of the associated CO 
clouds is well correlated with the Spitzer 8.0 and 24 jim emis- 
sion, and there are two giant elongated molecular clouds (AD 
and BC) of CO emission appearing on opposite sides of the ring- 
like shell of N131. The two elongated clouds may be triggered 
by the stellar wind from the clustered ionizing stars within the 
N131, but this possibility needs to be explored further. 

2. Contrasting the morphologic distributions between the CO 
and 8.0 /mi emissions, we found that the CO radii agree rather 
well with those from the 8 yum image. Dni3i = 8.6 kpc was 
adopted as the distance of the bubble N131 from the Sun. 

3. We found there is a cavity within the bubble at 1.4 GHz 
and 24 /mi, indicating that hot dust and ionized gas have likely 
been evacuated by stellar wind. 

4. Seven IRAS point sources (IR2, IR3, and IR8) are dis- 
tributed along the ringlike shell of the bubble N131. IR2 (IRAS 
19499+2613) is located at the peak of HII region G63.049- 
0.348. IR2, IR4, IR6, and IR8 are well correlated with the molec- 
ular clumps A, H, E, and G, respectively. 

5. We found 15 ionizing stars and 63 YSO candidates. The 
clustered YSOs are distributed along the elongated clouds AD 
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Fig. 6. Integrated intensity contours of the 12 CO emission of the blueshifted and redshifted clouds superimposed on the 1.4 GHz 
NVSS continuum emission contours. The contour levels range from 3.20 to 19.20 by 1.60 K for the blueshifted cloud, and from 
3.86 to 23.17 by 1.93 K for the redshifted cloud. The integration range is from -14.5 to -1 1.5 km s" 1 for blueshifted cloud, and from 
-9.5 to -6.5 km s" 1 for redshifted cloud. The contour levels of the 1.4 GHz NVSS continuum emission are 0.95, 1.26, 1.58, 3.16, 
6.32, 12.64, 25.29, and 50.58 mJy beam -1 . G63. 049-0. 348 is an HII region. The letters from A to I indicate the positions of nine 
molecular clumps, and the symbols "■" indicate the positions of IRAS point sources. The green ellipse indicates the position of 
bubble N131. The beam size of 12 CO emission is given in the filled black circle. 



and BC, and several YSOs are located around the clumps G and 
I. These locations may be the birth places of star formation trig- 
gered by the bubble Nl 3 1 . 
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Fig. 7. (a) GLIMPSE CC diagram [5.8]-[8.0] versus [3.6]-[4.5] for sources within a circle of about 10' in radius cente red on N131. 
The classification of class I, II, and III indicates different stellar evolutionary stages as defined bv lAUen et all (12004 . (b) 2MASS 
CC diagram (H - K s ) versus (/ - H). The sources for classes I and II are these, detected simultaneously by /, //, K s , 3.6, 4.5, 5.8, 
and 8.0 [i bands; the sources for class III are located within a circle of 2.60' (0.5xr^) in a r adius centered on N131. The gray squares 
represent the location of the main sequence and the giant stars (iBessell & Brett! [l988b . The parallel dotted lines are reddenin g 
vectors. The adopted interstellar reddening law is A j/A v = 0.282 , A H /A V = 0.175, and_A #, /A v = 0. 1 12 dRieke & Lebof skv 1985), 
and the intrinsic colors (H - K s )o and (/ - H)o are obtained from Martins & Plez (2006). 
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